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We present approximate solutions for the gluon and ghost propagators as well 
as the running coupling in Landau gauge Yang-Mills theories. We solve the 
corresponding Dyson-Schwinger equations in flat Euclidean space-time without 
any angular approximation. This supplements recently obtained results employ- 
ing a four-torus, i.e. a compact space-time manifold, as infrared regulator. We 
confirm previous findings deduced from an extrapolation with tori of different 
volumes: the gluon propagator is weakly vanishing in the infrared and the ghost 
<D propagator is highly singular. For non-vanishing momenta our propagators are 
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in remarkable agreement with recent lattice calculations. 
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The non-perturbative structure of the gluon propagator is of fundamental interest in 
Quantum Chromo Dynamics, and its knowledge would provide an important input for many 
calculations in hadron physics |L|. As infrared singularities might be encountered, a non- 
perturbative continuum method seems mandatory. In this context, the Landau gauge Dyson- 
Schwinger equations have been solved analytically in the infrared, but had to be truncated in 
order to get a closed system of equations |3|, f|, |5|, ffl. On the other hand, lattice calculations 
j7|, [| include all non-perturbative physics of Yang-Mills theories but are limited for small 
momenta by the finite lattice volume. Despite these shortcomings both approaches agree 
surprisingly well in their general statements: there is clear evidence for an infrared finite or 
even vanishing gluon propagator and a diverging ghost propagator. This is in accordance 
with the Kugo-Ojima confinement criterion, which in Landau gauge includes the statement 
that the ghost propagator should be more singular than a simple pole 0. 

One of the obstacles encountered in extracting definite values of infrared exponents in 
Dyson-Schwinger studies is the angular integrals inherent to these equations. Therefore 
approximated treatments of the angular integrals have been applied in numerical works so 
far [0, [3] . In general the angular approximations proved to be good for high momenta but less 
trustworthy in the infrared. Recent studies |4|, [5], |6| therefore concentrated on the infrared 
analysis, where exact results have been gained for the limit of vanishing momentum. 

In a previous paper [ ID | a truncation scheme has been presented, which on the one hand 
provides the correct anomalous dimensions of the ghost and gluon dressing functions, Z(k 2 ) 
and G(k 2 ), in the ultraviolet region of momentum. On the other hand, it reproduces the 
infrared exponents found in || § which are close to the ones extacted from lattice calculations 
[0, ||. The coupled set of equations for the ghost and gluon dressing functions on a four- 
torus have been solved recently flO |. This turned out to be an effective tool to overcome 
the angular approximation for non-vanishing momenta. The numerical solutions proved to 
be compatible with only one out of two solutions of the infrared analysis of refs. |5|, 
thus suggesting that not every analytical solution for zero momentum connects to numerical 
solutions for non- vanishing momenta. However, employing the torus as an infrared regulator 
shares with lattice calculations the problem of working on a finite volume. Therfore further 
investigations are highly desirable. 



In this letter we present the solution of the coupled set of equations described in ref. (10| 
in flat Euclidean space-time. The results confirm our previous findings on the torus. We 
obtain a numerical solution for the gluon and the ghost propagator with infrared behavior 
D g i(k 2 ) ~ {k 2 ) 2K ~ l , k m 0.595 and D g h(k 2 ) ~ {k 2 )~ K ~ l . Correspondingly we find an infrared 



1 



.1 .1 



C7> 



Figure 1: Diagrammatic representation of the truncated gluon and ghost Dyson-Schwinger 
equations studied in this letter. Terms with four-gluon vertices have been dismissed. 

fixed point of the running coupling, a(0) ~ 2.972 for the gauge group SU(3). 

The truncation scheme 

To make this letter self-contained we will give a short summary of the truncation scheme 
introduced in detail in ref. HlQ| . It has turned out in previous studies 0, that the dressing 



of the three-point vertex functions does not alter the qualitative behaviour of the ghost and 
gluon system. We therefore restrict ourselves to bare vertices for simplicity. Furthermore all 
diagrams including four-gluon vertices have been dropped. A diagrammatical representation 
of the resulting system of equations is presented in Fig. [I]. To extract the gluon dressing func- 
tion the gluon equation has to be contracted with a tensor. Maintaining Lorentz invariance 
its general form reads 

rjS(k) = 8 lw -c^. (i) 

The parameter £ allows to continuously interpolate from the transversal projector, ( = 1, 
to the Brown- Pennington projector JTTJ] , ( = 4. In a full treatment without truncation the 
gluon equation would be transverse and the propagators would not depend on (. We will 
see later on to what extent transversality is violated in the present truncation scheme. 

Closely connected to the subject of transversality is the appearance of unphysical qua- 
dratic divergencies in the gluon equation. For ( = 4 the equation is projected on the longitu- 
dinal tensor k^k u /k 2 and no quadratic ultraviolet divergences appear |lTJ. For a transverse 
polarization the gluon equation would be independent of (, and thus no quadratic diver- 
gences would occur for all ( in a full equation. In the practical case, the loss of transversality 
causes quadratic divergences to appear and their careful removal is essential for a numerical 
solution of the system. However, this procedure should not change the infrared properties 
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of the equations. Anticipating that the ghost loop is the leading contribution to the infrared 
we subtract the quadratically ultraviolet divergent constant in the gluon loop. 

The resulting system of equations for the ghost dressing function G{x) and the gluon 
dressing function Z(x) reads 

1 - 9 9 2 N C [ ^^(y)Z W , (2) 



G(x) ° * J (2tt) 4 xy 

1 7 , 2 N c [ d A q M(x,y,z) 2 N C f d A q Q(x,y,z) 

W) = 3 + 9 t J Wf^~ G{y)G{z) + 19 t J - z{v)z{z) ■ 

(3) 

The symbols x := k 2 and y := q 2 denote squared momenta where z = x + y — 2^fxy cosO. 
The kernels ordered with respect to powers of z have the form: 

(x ~y) 2 \ If x + y \ _ 1 
z 2 V 4 J z V 2 I 4 



/\ (./•. //. : j = ±(-S±^)+t(±^)-t j (4) 



M(x,y,z) = - £ H H , 5 

v ' y ' ; 2 V 4 2 4 x J 2 2x 4x K ' 

1 A^ 3 2 19 ~C 5-C 2 Cz/ 3 
Q{x,y,z) = —1-—+X — xy + —— y + - — 

z l \ 8 y 8 4 8 i 



i fx 2 i5+c i7-c M y 

+- — x — y + C— 

z \ y 4 4 x 

19-C s , 17-C , 9CZ/ 
8 y 4 Ax 

+ <§ + 5 ^) + ' s 4 + J (4 - - (6) 

Here the last term in the gluon kernel Q is introduced by hand to subtract all quadratic 
divergencies from the ghost and gluon loop, see ref. PD[ . 

In order to ensure the correct one loop scaling of the dressing functions in the ultraviolet 
within this truncation scheme one has to violate the Slavnov-Taylor identity Z\ = Z3/Z3, 
similar to the treatment in |2| . Refering the reader to the detailed derivation in |7J , we only 
give the resulting momentum dependent substitution Z\ for Z\\ 

Giy) (l-a/S-2 a) G{z) (l- b/ S-2 b) 

Z 1 ^Z 1 (x,y,z)= z{y){1+a) z{z){l+b) (7) 

With this substitution the logarithmic behaviour of the gluon loop in the ultraviolet is the 
same as the one of the ghost loop. The coefficients in the kernels then sum up to the 
correct anomalous dimensions, namely 5 = —9/44 for the ghost and 7 = —13/22 for the 
gluon propagator. The parameters a and b can be chosen freely. In our calculations we use 



b = 35, which has been shown |10| to minimize the momentum dependence of Z\. Note 
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that a = b = leads to the truncation scheme of ref. [12| and a = 35, 6 = together with 
appropriate vertex dressings to the one of ref. @. 



Infrared analysis 

In the truncation schemes of refs. [|, |5], ^, the ghost loop provides the infrared 



leading term in the gluon equation. The only exception is given in ref. ||12|| : in this scheme 
also the gluon loop contributes to the same order as the ghost loop. As we do not want to 
change the infrared behaviour by manipulations of the renormalization we choose a = b = 35 
which is the only choice to give a constant Z\ in the infrared. Then the gluon loop does 
not contribute to the infrared behaviour of the equations. Therefore the leading analytical 
solutions found in the ghost loop only truncations of [[| with £ = 4 and of |5], |6| with ( = 1, 
are found in our case as well [|Kj] . In the following we will determine subleading contributions 
in the infrared, and judge their importance for the numerical treatment of the coupled system 
of equations. 

We employ the ansatz 

Z{x) = Ax 2K (l + fx p ) 

G{x) = Bx- K (l + gx p ) (8) 

in the eqs. (0) and (^). The left hand side and the integrands are expanded to appropri- 
ate powers of momentum. After integration (for technical details see refs. || |10, |l3j) the 



coefficients of equal powers on both sides of the equations have to match for consistency. 
The conditions on the leading term remain unchanged. Hence for an arbitary projector the 
leading power k is given by 

(2 + k)(1 + k) _ 4k -2 = 

18(3 -2k) 4Ck - 4k + 6 - 3C 

This equation has at least two solutions for a given £ from which only one connects to the 
numerical solution for finite momenta ||10|| . The result for the running coupling is given by 

a(0) = ^ r(3 - 2/t)r(3 + /t)r(1 + K) ^ 2 972 (10) 
a[()) 18 ^(2 - k)T(2k) ' 1 j 

where the nonperturbative definition a(x) = a(fi 2 ) Z (x)G 2 (x) has been used. 

Matching subleading powers leads to the coupled set of homogeneous equations for /, g 
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and p: 

6«(k-1)(-3 + 2k) r(2«)r(2-« + p)r(3-«-p) 

'« + p - 2)(« + p - 1)(« + p) r(4 - 2«)r(i + « - p)r(3 + « + p) ~ 1 ' u 

f z 3(-2 + 2k + p) r(2 - K )r(2K + p + l)T(3-K-p) \ 

V Z/ 2( K + p - 2)(k + p - i)(k + p) r(3 - 2k - p)r(i + «)r(3 + « + p) ; 7 

4CK-4K + 2p-2Cp + 6-3a r(2 - K )r(2 K - p)r(2 - K + p) \ 

V 2K-p-l r{l + K)T{l + K-p)T(i-2K + p)J 9 1 ' 1 ' 

Here v = N c g 2 AB 2 /48ir 2 = ot(0)/47r. There is either the trivial solution / = g = or one 
has to set the determinant of these linear equations to zero. We then obtain the results 
P(i) = 0, p(2) = 0.58377, p( 3 ) = 1.20300 and several other solutions with higher values of 
the power. The solution p^ = corresponds to the pure power solution. The lowest 
nonvanishing solution, p^) — 0.58377, is sufficently high that we safely may neglect it in the 
numerical treatment of the infrared part of the equations. This will be detailed in the next 
section. 

Numerical Method and Results 

Two of the four integrals in equations @ and (^|) can be done trivially and yield a factor 
of 47r, while the remaining angular integral and the radial integral have to be performed with 
the help of numerical routines. To achieve high accuracy we split the radial loop integral 
into three parts, y G [0, e 2 ], y G (e 2 ,x] and y G (x,xjjv]- The integrals will be computed 
numerically, however according to the value of their argument the dressing functions Z and 
G have to be handled differently. In the infrared region, y,z G [0,e 2 ], Z and G behave like 
powers and are replaced according to equations (H) with f — g — 0. The infrared matching 
point is chosen to be e 2 = (0.55MeV) 2 in our calculations. In the high momentum regime, 
y G (x,xuv), arguments z occur which are larger than the numerical cutoff xuv- There we 
approximate the respective dressing functions by the expressions 



Z z) = Z s) 



oo\og (~) +1 
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G(z) = G(s) [wlog(^) • (13) 

according to the one loop behaviour of the solutions as has been detailed in [|Kj. Here to = 
llA^ c a(s)/127r and the squared momentum s is a perturbative scale. We chose s = (174GeV) 2 
to be slightly lower than the numerical cutoff xjjv = (177GeV) 2 . To be able to perform the 
angular integrations for momenta [e 2 ,X[/y] we expand the dressing functions in Chebychev 
polynomials and solve the coupled system of equations for the expansion coefficients using 
a Newton iteration method. Details of this technique can be found in ref. H. 
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Figure 2: Shown are the results for the gluon dressing function, the ghost dressing function 
and the running coupling using a transverse projector. The solutions are compared to pre- 
vious results on a torus, for the running coupling two results for two different volumes are 
presented. The dashed lines are the fits given in eq. (|14|). 



We apply a MOM regularization scheme similar to the ones used previously in refs. 0, [3] . 
For numerical reasons it is favourable to subtract the ghost equation at zero momentum and 
the gluon equation at the scale s introduced above. Due to this the unknown renormalization 
constants Z 3 and Z 3 drop out and instead of them one has to specify two input variables. We 
chose the infrared parameter A (c.f. eq. (0)), and implement the condition Z(/z 2 )G 2 ( / u 2 ) = 1. 
Furthermore, one has to specify the coupling a(fi 2 ) = g 2 /lQTr 2 entering eqs. (0,0). Finally 
we have fixed the momentum scale by calculating the running coupling for the colour group 
SU(3), and requiring the experimental value a(x) = 0.118 at x = M\ = (91.187GeV) 2 @. 

In Fig. we show the results for the transverse projector, ( = 1. The gluon and ghost 
dressing functions behave powerlike for low momenta and obey one loop scaling in the 
ultraviolet as expected. Accordingly the running coupling has a fixed point in the infared 
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Figure 3: Shown are the results for the gluon dressing function, the ghost dressing function 
and the running coupling, c.f. fig. 0, for different projectors. 



and decreases logarithmically in the perturbative regime above several GeV. For intermediate 
momenta the behaviour of the running coupling seems to suggest a second zero of the (3- 
function around (lOOMeV) 2 . However, such an extremum in a(x) would result in a double 
valued /5-function, therefore we regard it as an artefact of our truncation scheme. In fact 
with subleading gluon contributions in the infrared and perturbative suppression of the two 
loop contributions in the ultraviolet it is exactly the intermediate momentum regime where 
the omission of the two loop diagrams should cause deficiencies in our solutions. 

The asymptotic behaviour of the solutions can also be seen from the functional form of 
our fits. The dressing functions Z(x) and G(x) can be described by 

a(0) „, x cx K + dx 2n 



a(x) 



Z(x) 



ln(e + aix a2 + bix b2 ) ' 
I \ \ 1 + 25 



R(x) 
G(x) 



1 + cx R + dx 2n ' 

8 

R-Hx) 



a(x) 



(14) 
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with fitting parameters a 1; a 2 , b\, b<i and c, d for the running coupling a(x) and the auxiliary 
function R(x) respectively. The values of a(0) = 2.972 and k = 0.5953 are known from the 
infrared analysis in the previous paragraph, whereas 5 = —9/44 is the anomalous dimension 
of the ghost dressing function and a(fi 2 ) = 0.9676. The six parameters of the fit are given 
by a x = 5.292GeV" 2a2 , a 2 = 2.324, b t = 0.034GeV" 262 , b 2 = 3.169, c = 1.8934GeV" 2K and 
d = 4.6944GeV _4K . As can be seen in Fig. |2] the fit works very well and can be used as input 
for phenomenological calculations in future work. 

Our results for different values of the parameter ( (c.f. eq. (|1])) are shown in Fig. |3|. 
In accordance with the infrared analysis the power k changes from k = 0.5953 for £ = 1 
to k — 0.4610 for ( = 5. The perturbative properties of the solutions remain unchanged. 
The bump in the running coupling gets smaller but does not disappear even for ( = 5. It 
has already been stated above that the dressing functions would be independent of £ in a 
complete treatment of the gluon equation. As all our solutions are very similar even on a 
quantitative level we conclude that transversality is lost only to a moderate extent. This is 
a somewhat surprising result in such a simple truncation scheme as the one at hand. 

The Brown-Pennington projector, ( = 4, is an exceptional case as can be seen from 
eq. (|9|). Here the ^-dependence of the second term cancels and only one solution, k = 1, can 
be found {c.f. ref. H). We found no numerical solutions for this case. However, within the 
limit of numerical accuracy, solutions for ( slightly different from 4 can be found leading to a 
value for k slightly different from 1/2. E.g. in Fig. |^ the case ( = 3.9 leading to k = 0.5038 
is depicted. 



Finally, we compare our results to recent lattice calculations for two colours [15|. As 
the solutions on the lattice include all non-perturbative effects, the results shown in Fig. |] 
confirm again that the omission of the two loop diagrams in our truncation scheme mostly 
effects the region around the bending point at 1 GeV. However, in general the qualitative 
and partly even quantitative agreement of the two methods is remarkable. The combined 
evidence of the two methods points strongly towards an infrared vanishing gluon propagator 
and an infrared singular ghost propagator in Landau gauge. 

Conclusions 

In this letter we have presented approximate non-perturbative solutions for the gluon and 
the ghost propagators as well as the running coupling in Landau gauge. We obtained these 



solutions for the Dyson-Schwinger equations in the truncation of ref. |10|, working with 



bare vertices and omitting all diagrams that involve four-gluon vertices. An improvement 
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Figure 4: Solutions of the Dyson-Schwinger equations compared to recent lattice results for 
two colours [131. 



to previous treatments has been the explicit numerical calculation of all angular integrals 
thus overcoming the angular approximations that have been made so far. It is interesting to 
note that only one solution (for a given projector) has been found. Furthermore, we could 
demonstrate that our truncation scheme violates transversality of the gluon equation only 
to a moderate extent. Previous findings on a four-torus have been checked and found to 
be in agreement with the flat space-time solutions up to effects of finite volume in the very 
infrared. Despite the simplicity of the truncation our solutions agree remarkably well with 
recent lattice calculations. For the infrared behaviour of Landau gauge Green's functions 
circumstancial evidence points towards a weakly infrared vanishing or even finite gluon 
propagator, while the ghost propagator seems to rise more strongly than a simple pole. 
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